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ON SAMPLING OF STATIONARY INCREMENT PROCESSES 
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Under a complex technical condition, similar to such used in ex¬ 
treme value theory, we find the rate g(e) -1 at which a stochastic pro¬ 
cess with stationary increments £ should be sampled, for the sampled 
process £(L'/9( e )J ?( e )) to deviate from £ by at most e, with a given 
probability, asymptotically as e J. 0. The canonical application is to 
discretization errors in computer simulation of stochastic processes. 


1. Introduction. Let {£(i)}tg[o,i] be a stationary increment (SI) process. 
In Section 4, under Condition © from Section 2, we prove that, for suitable 
functions q\, q 2 , w, F > 0, a constant k > 0 and an open interval J C M, 
denoting = sup{/c G Z : k < x}, 


( 1 . 1 ) 


( 1 . 2 ) 


lim p r supteio,!]^^) — g(L^/gi(g)Jgi(g))] -g 

e|0 \ w { £ ) 


limpf SUPte[0 ’ 1] ^ — ^(L*/g2 (g) Jg2 (g))| -g 

ej.0 \ w ( £ ) 



for x G J. 

= g -UF(x) 


for x G J. 


Methods of proof relate to Albin (1998), for the analysis of single sam¬ 
pling intervals, and to Albin (1990), when assembling these to the limits 
(1.1) and (1.2). The main inspirations for those works in turn were Berman 
(1982, 1992) and Pickands (1969a, b) together with Leadbetter, Lindgren 
and Rootzen [(1983), Chapters 12 and 13], respectively. 

In Section 3 we discuss Condition © and the role of (1.1) and (1.2) in 
applications. 
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Versions of (1.1) and (1.2) have a history for Gaussian processes; see 
Belyaev and Simonyan (1979), Piterbarg and Seleznjev (1994), Seleznjev 
(1996) and Hiisler (1999). Section 5 features applications to stable processes. 
Future applications could be to light-tailed Levy processes and Ornstein- 
Uhlenbeck processes [cf. Albin (2003a, b)]. 

Albeit similar in appearance, (1.2) has little if any relation to the small 
ball problem, to estimate P{sup tG r 0 n |£(i)| < e} as e J, 0 [e.g., Kuelbs and 
Li (1993)], as (1.1) is unrelated to estimates of P{sup tg r 0jl i £(i) < e} [e.g., 
Marcus (2000)]. 


2. Condition ©. The conditions (2.1)-(2.8) on the process {^(i)} te [o,i] 
listed below will be collectively referred to as Condition ©(*), with i = 1 or 
i = 2: 


(a) Let q = q(e),q = q(s),q = q(s) and w = w(e) be strictly positive func¬ 
tions on (0,oo), with q{e),w{e) j 0 as e j 0, and with either q = 1 or Q\ = 
lim e j_o 1/? = 00 , and either q = q or Q 2 = lim e j_o q/q = 00 . Set I = (0, Q 1 ) or 
7 = (0,Qi]. 

(b) Let £ be SI, separable and measurable, defined on a complete prob¬ 
ability space. Depending on whether (1.1) or (1.2) is in focus, put £ e = £i )£ 
or ie = 6 )£ , where 


£i, £ (*) = [£(<#) -£0?I_£-OJ) - e]/w 

= [|^(9*)-^(?L*-°J)I -e]/w 


for t > 0. 


(c) Pick an open set J Cl with 0 E J, a continuous function F : J —> 
(0, 00 ) with F{ 0) = 1, and an integrable function G : [0, Q\) —> (0, 00 ) with 
G(0) = 1. Assume that, for 0 < S/Q 2 < Qi — r < Q\ and x E J, 


( 2 . 1 ) 


lim 

ej.0 


P{£ £ (1 — rq — sq ) > x} 


(■ q/q)q 

(d) Assume that, for T G I and x G J, 

TA(i/g) P{£ e (l — rq) > x} 


= F(x)G ( r + — 


( 2 . 2 ) 


lim 

eJ.o Jo 


Pfe(l)>0} 


dr 


= F(x) / G(r) dr = F(x)E(T ; G) > 0. 


(e) Assume that, for n GN, r G (0,Qi), t £ ( 0,rQ2) n and (independently 
of) x £ J, 


( 2 -3) l £ ifo P |n{^(l -rq + Uq) > x} £ e (l -rq) > x | =P r (t). 
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(f) If Q 2 = oo, assume that, for any x G J, a.e. for r£/, 

rrq/q 

(2.4) Jim limsup / P{£ e (l — qr + gs) > x|£ e (l — qr) > x} ds = 0. 

T—>oo efO JT 

(g) Assume that, for T E / and in J, with f denoting the direct 

upper Riemann integral of nonnegative functions [e.g., Asmussen (1987)], 


1 


(2.5) 


lim lim sup ■ „ r _. 

40 ej.0 «P{?e(l)>0} 

fTA(l/q) 


<p Ps sup £ e ((l - tq-saq) + ) > x, 

Jo Ue[o,i] 

Ce(l - tq) V 4((1 - tq- aq ) + ) < x j dt = 0. 


(h) If I = (0, Qi), assume that, for Tel and x G J, 

( 2 . 6 ) 


• f limsu _ P {^Pr S [TA(i/ g -),(i/ g -))^(l- rg)>x} _ 

S7 m) P{e £ (i)>0} 


(i) Assume that, for T £ 7 and x G J, 


LVgJ 


(2.7) lim lim sup sup ^ 


?P{^(1 - r<?) > *;4(ra - sq) > x } 


40 £i o r,*e[0,TA(i/9)) ^ PfeC 1 “ r{ ?) > 0}P{^(1 - sg) > 0} 


= 0 . 


(j) Given T € I , x E J and a, h > 0, assume that, for any s\,..., Si,t\,..., t# 
[0,1/g] n(J^ =1 (n —TJ, n\ with Sj + \ — Sj > aq, tj + \ —tj > aq and t\ — Si>h/q, 


p \ D {^(s k ) < x}. p| {te(tt) < x} 


ej.0 


( 2 . 8 ) 


.fc=1 


£=1 


-P< p{^(sfe)<a;}>P< P {€e(te) < x} \ =0. 


,fc=i 


i=i 


Theorem 1. If the Condition &{i) holds, for an i E {1,2} ; then (l.i) 
holds for some constant k > 0. 


3. Discussion. The idea behind (1.1) and (1.2) is to use extreme value 
technology, for the extreme values {s + xw} xe j of the processes £(qt) — 
£,(q\t — 0J) and | £(qt) — £(g[i — 0j)|, where q = q{e) is chosen small enough 
to make e an “extremely large” value for these processes as e j 0. Conditions 
(2.1)-(2.8) emerge from inserting £(qt) — £(q[t — 0J) and | £{qt) — f,(q[t — 
0])| in analogous conditions from extremes [e.g., Albin (1990, 1998)]. What 
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makes life difficult here is that these processes depend on the level e and do 
not possess any invariances like self-similarity (SS) or SI. 

However, for £ SSSI, (1.1) and (1.2) are easily transformed to more con¬ 
ventional extreme value problems, involving processes that do not depend 
on the level. 

Schemes from extremes are natural take-off points to check Condition (5. 
Basically, we refer readers to the literature on extremes for that purpose (but 
see Section 5). In fact, very much can be written on how to check (2.1)-(2.8) 
[see, e.g., Albin (1998) for a hint on this]. We plan to return to this in the 
future. 

We discuss Condition &(i) for i = 1; the case when i = 2 is very similar. 

Condition (2.1) implies that 


(3.1) P{&(l)>0} = P{£( g )-£(0)>e}~(g/g) g as e 10 


and that 



(3.2) 


F(x). 


This means that £(q) — £(0) belongs to a domain of attraction of extremes 
[e.g., Albin (1990), Section 1], modulo the fact that this random variable 
depends on e. In addition to (3.1) and (3.2), (2.1) requires that q is chosen 
in such a way that 



where GeL 1 . 


Arguably, in the presence of (3.2), this is a reasonable requirement. 

In the presence of (2.1), condition (2.2) is just a technicality. Notice that, 
by (2.1) and Fatou’s lemma, (2.2) for a certain T = T implies (2.2) for T <T. 
Typically, (2.3) is verified by means of establishing that 


(3-4) {(€e(l-rq + tq)\£ e (l-rq) > z)} t > 0 4 {Cr(t)}t>o as e | 0, 


for some stochastic process {Cx(^)}t>o- This together with (2.1) implies (2.3), 
since 


limsupP{£ e (l — rq + tq) = x|£ e (l — rq) > x} 

ej.0 




In particular, (2.3) implies that 


(3.5) limP{£ e (l - rq + tq ) > x|£ e (l - rq) > x) = P r (t). 

eJ.0 
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Assuming that (2.3) is established via (3.4), a formal argument gives 

rrq/q 

lim / P{£ e (l — qr + qs ) > x | £ e (l — qr) > x} ds 
eio Jt 

poo r poo 'j 

= P{Cx(s)>x}ds = E^J^ l( 3; ,oo)(Ca ; (s))dsj 

when Q 2 = 00 : the interpretation of (2.4) is thus one about asymptotically 
small mean sojourn times above the level x, for the process C, x at infinity. 

Condition (2.5) is typically verified by bounds on increments of £ e (see, 
e.g., Section 5). For example, (2.5) holds if, for some constants K, d > 0 and 
e > 1, 

(3.6) P{£ e (l -rq + sq) < x -5, £ e (l - rq) > x} < K5~ d s e P{£ e (l) > 0} 


for e large enough, with uniformity in r. When (2.1) and (3.4) holds, we 
have 

V P{Ce(l -rq + sq) <x -5, £ e (l - rq) > x} 

l-m- P{fe(I)>0 } -= F WG(r)P{CxW <x-S). 

For F and G “nice,” the interpretation of (3.6) thus is [recall that Cr(0) > x] 

P{C x{s)<x-5}<K5- d s e . 


Condition (2.5) is often easier to verify for I = (0, Q\) than for I = (0, Qi]. 
The price for this is that (2.6) is required, to take care of £ e (i) for t close to 
zero. 

Condition (2.7) is a mixing conditon. It requires in essence that 
P{£e(l - rq) > x,£ e (n - sq) > x} < A"P{£ £ (1 - rq) > 0}P{£ e (n - sq) > 0} 


for some constant K > 0. Condition (2.8) is also a mixing conditon that is 
imposed for process values that are farther apart than those in (2.7). 

In most applications, Q 2 = 1 so that q = q, since Q 2 = 00 usually corre¬ 
sponds to rather pathological processes: for example, for fractional Brow¬ 
nian motion, one has Q 2 = 00 for the index of self-similarity H E (0, ^), 
while Q 2 = 1 for H E [^, 1) [see Piterbarg and Seleznjev (1994) and Selezn- 
jev (1996)]. 

For Q 2 = 1, the functions q, w and q = q are typically determined from 
(3.1), (3.2) and (3.3), respectively (in that order). For Q 2 = 00 four functions 
are sought, and an additional equation is required, for example, (3.5). Of 
course, a rigorous use of Theorem 1 requires a verification of the entire 
Condition ©. 

The following two examples are canonical for the role of (1.1) in applica¬ 
tions: 
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Example 1. In a computer simulation a largest deviation d = e + xw{e) 
between the discretized simulated process £(L-/<z(£)J <?(£)) and the process 
{£(*)}te[o,i] itself is tolerated, with a certain probability p = e kF ( x ') . For this, 
the sampling rate should be g(e)” 1 , asymptotically as e j 0 (same thing as 
d | 0). Here the deviation is measured “two-sided” in (1.2), and “one-sided” 
in (1.1). 

An applied risk application can be to estimate a high quantile 

(3.7) P< sup £(t) > u > < 2p for u large and p small. 

Ue[o,i] J 

This is solved using computer simulations, to find a u such that 

rUM 6 )! 1 _ 

P< [J {£(kq) > u — d} > = p = 1 — e~ nF ^ where d = e + xw(e), 

l k= 0 J 

for e small enough that the limit (1.1) is accurate, because then (1.1) gives 


P 


sup £(t) > u 

*€[ 0 , 1 ] 


< P 



u — d} 


+ P 


sup [£(t) -^(lt/g(£)Jg(e))] >e + xw(e) 

*€[ 0 , 1 ] 


< p + 1 - e“ KF(3:) = 2 p. 


Example 2. In a theoretical risk application, the problem is again to 
solve (3.7). Let £ be stationary, with £(0) belonging to a domain of attrac¬ 
tion, that is, 


(3.8) 


P{£(0) >u + yw(u)} 
P{£(0 )>u} 


= H(y) 


for y € J, 


for some functions w,H > 0. Picking a — y E Jn (—oo,0), (3.7) is solved 
taking 

yw(u) = e + xw(e) and + l^i7(—y)P{£(0) > u} =p = 1 — e~ KF ^ x \ 

for e small enough and u large enough that the limits (1.1) and (3.8) are accu¬ 
rate, because then we have, for such e and u, by stationarity, (1.1) and (3.8), 


P< sup £(t) > u 
Ue[o,i] 
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( U/9( £ )J 

< p \ U {£(kq) > u - yw(u)} 
l k =0 


+ P 


sup 

te[o,i] 


m - c 




9(e) 


> e + xw{e ) 


< + l) P{£(0) >u- yw(u)} + 1 - e~ nF ^ = 2 p. 

We conclude this section with an example featuring, possibly, the only 
nontrivial process for which the probability in (1.1) can be calculated ex¬ 
plicitly. 


Example 3. Let B be Brownian motion (arbitrarily started), W = B — 
B{ 0) and 


9*00 


; 1 ■ i 


+ In 4 In 


| 21n / 2Var{W(l)}i 

V V27T 


1 -1 


Var{W(l)} 

By self-similarity and other elementary properties of Brownian motion, we 
get 

t 


P<{ sup B{t) — B 

[o,i] 


L <7 J 


9 <e 


= ( P<j sup W(t) < e 
^e[o,g] 


X) ‘ p{ sup W(t)<e\ 
J) Ue[o,i-Li/gJ q] J 

(1 - 2P {W(q) > £}) [1/qi (l - 2P| w(l 


1 

L9J 


q >£ 


, y/2C~q 

1--=— exp 


2 Cq 


(l + o(l)) 


1/9 


-2 


4. Proof of Theorem 1. The sojourn time 

rsA(l/q) 

L e (s;x) = J l(x )0o) (Ce(l -qr))dr, ( s,x)elxj , 

satisfies 

(4.1) E{L £ (T;x)}~P{e £ (l)>0}P(x)E(T;G) as e 10 

for T E I and x E J [by (2.2) and Fubini’s theorem]. By the elementary 
identity 

j'T/\(l/q) 

l(y,oo) (■^e(-^"'> ■1')) / 1(— oo,y] (-^e(®) -e)) 1(x,oo) (£e(l qs~)) ds 

J 0 

l(y,oo) %)^)y 



















J .M. P. ALBIN 


for T G I and y > 0, it therefore follows that 
If P^qL e (T;x)/q> z^dz 

roo 

= P {(L £ (T;x)-qy/q)> z}dz 

Jo 

(4.2) =E{(L £ (T-x) - qy/q)ltf y /q !00) (L e (T;x))} 


E'j ^~{qy/q, oo)(E £ (T, x')') 

r TA(l/g) 


fTA(l/q) _ 'j 

J [1 1(— oo,qy/q] (-^e (®> ®))] ^-(x,oo) (£e (1 (7®)) J 


= E 


rT A(l/?) 


l(qj//g,o°) 2'))l(x,oo) (^(1 (js)') ds 

rTA(l/q) 

/ P{^e(s;*)/9>y>Ce(l-9s) >x}rfs. 

JO 

Using (2.1) and (2.2), we get density convergence 

for s G [0,T) 


lim P{^(1 - qs) > x} _ G(s) 


U0 E{L £ (T; x)} E(T;G) 

Hence (4.2) together with Scheffe’s theorem [Scheffe (1947)] show that, as 

£|0, 

p {qL e (T-,x) > qz} 


f 

J V 


■ dz 


(4.3) 


(q/q)E{L £ (T ]X )} 

( r s q/q 1 

y o l( 3 ;,oo)(^(l -qs + qt )) dt > y £ e (l - qs) > xj 

P{4(1 ~qs)>x} 


r T A(l/g) ( f s q/q 

p< 


■ ds 


L 


TA(l/q) ( fsq/q 

P 


E {L e {T-x)} 

f fsq/q 'j 

l(i,oo) (£e(1 -qs + qt))dt >y £ e (l ~ qs) >xj 

G(s ) 


ds. 


E{T-G) 

Picking a T > 0, (2.3) gives convergence of moments 

r{sq/q)AT 

/0 


r / r(sq/q)AT \n >, 

E| (jf 1 (z,oo)(&( 1-<7'S + ^))^J 6(l-^)>xj 
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'[0,(sg/g)AT]' 


nw-v + to) ^ ~Q s )> x \dt 


< 1=1 


•'[0,(sQ 2 )AT]" 

Here we have, with obvious notation, 

r(sq/q)AT 


P s (t ) dt as £ | 0 for n E N. 


/ rtsq/q)A1 \ 

(j l(a:,oo)fe(l - < 7 s + gt))df £ e (l -gs) >xj 


—> some random variable C ~ 

^S,I 

(since the left-hand side is bounded by T). Using this in (4.3), we get 

1 P {qL e (T;x) > qz} 


lim inf 


(4.4) 


UO J y (q/q)E{L e (T-,x)} 
> 


dz 


TA(l/q) Q( s ) 

P K,t >y } e(T-G) ds 


= Af( r ;y) T A(T;y) asT|oo. 


On the other hand, using (2.4) if Q 2 = 00 , by (4.3) and Markov’s inequality, 
P {QL e (T;x) > qz} 


h “o UP 4 (q/q)E{L e (T-x)} 


dz 


< lim sup 
ejO JO 


rl A(l/ q) ( rsq/q 

/ P / ^l(x,oo ~ QS + qt)) dt > 5 

JO U(sq/q)AT 

c £ (l-gs)>x| 




E(T-G) 

r TA (V9) . . G(s) 

’e{T-,G) 


r-l /M.1/ q) 

(4.5) + lim sup / P{C^>y-d} 

ej.0 JO 

f T A(l/q) rsq/q P{£ £ (1 — qs + qt) > x|£ e (l — qs) > x] 
< Inn sup / / 

elO Jo J(s 


dt 


(sq/q)AT 


x^L-ds 

E(T-G) 


+ A dT-y-6) 


0 + A(T; y — 0) as T —> 00 and <5 j 0 (in that order). 
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By (4.5) and elementary arguments, we have, for T G I and x £ J, 


liminf ^ UPte[0 ’ TA(1 / , ‘)] ^ (1 ” ^ > X} 


ejO 


(4.6) 


1 

> - 


(q/q)E{L e (T ]X )} 

1 P {qL e (T-,x) > qz} 


1 — lim sup 


V L 

> lim sup 
yio 


10 Jy (?/?)E {L £ {T]x)} 
» sl » 


dz 


in a suitable way. Here the left limit is strictly positive, since by a version 
of (4.5), 

lim Slip T P{jL l (r i X L > 0 <fe 


10 Jy+x ( q/q)E{L e (T;x )} 

pTA(1/q) f p{sq/q)hT 

< lim sup / 
elO 40 


l(:r,oo )(6(l-^ + #))dt>y 

6(1-<?s) >®| 


G(s) J 
x — r ' ds 
E(T-G) 

r TA(l/q) rsq/q 


rl/\(i/q) rsq/q 

+ limsup/ / _P{6(1 — gs + qt) > x|6(l — qs) > x} dt 

elO JO J(sa/a)AT 


l(sq/q)/\T 


x-S£) * 

e(t,g) 

< 0 + - for T E (0, y] large enough. 

On the other hand, we have, for x < x < x in J and y > 0, by (4.4), as 
e| 0, 

p {sup te[ o,TA(i/g)] 6(1 - tq) > x} 

(' q/m{Le(T-,x )} 

1 


< 


y(<7/9)E{L £ (T;x)} 


rv 

x / P 


sup 6(1 - tq) > x \ U {qL e (T;x) > qz} 

0 ( Ue[0,TA(l/g)] J 

L((Tg)Al)/(a?)J 


LU4 N 

U (J {6(1 “ kaq) > x] > dz 

k= 0 J 













SAMPLING OF STOCHASTIC PROCESSES 


11 


1 

< - 


1 - 


f°° P{qL e (T-x ) > qz} 

y l Jy {q/q)V{Le(T-x)} 

L((T$)Al)/(a$)J 

+ E 


dz 


P {qLe(T;x) < qy,£e{ 1 - &ag) > 5} 


k=0 


+ 


1 


(9/g)E{L e (T;x)}- 


(9/g)E{L e (T;x)} 


sup £ e (l — tq)>x, 
^G[0,TA(l/g)] 


(4.7) 


Al)/(ag)J 'j 

H (£e(l ~kaq)<x}\ 
t=o J 


L((Tg)Al)/(ag)J 

r 

k= 


1 — A(T; y — 0) 

<--- + o(l) 


y 


+ 


p {qL £ (T-x) < qy,ie{ 1) > x} 
(q/q)E{L e (T-,x)} 


L((Tg)Al)/(ag)J fc 

+ E / [P{y L e( T ^) < &/, 

-to- 1 

£e(l - tog) v £ e ((l - (i + l)ag) + ) > x} 
x^/flE^T;®)})" 1 


+ P{6(1 - kaq) > x, 

£e(l - taq) V £ e ((l - (t + l)aq) + ) < x} 

x ((g/g)E{L £ (T;x)}) _1 ]dt 

1 

+ (g/g)E{L e (T;x)} 

L(( T 9)/\l)/(a l j)J 

x ^ P< sup £ e (l — taq)>x, 
k= o lte[fc,fc+i] 


£e(l - kaq) V £ e ((l - (fc + l)ag) + ) < x j 


< 1 - A(T ^~ 0) +o(l) 

y 
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+ lim sup 
40 


P {qL £ (T ; x) < qy , 6(1) > x} 


+ lim sup 2 

40 


0 q/q)E{L e (T-,x )} 

■LTA(i/g)j p{gL £ (i;x) < qy,£ e ( 1 - t<?) > x} 


0 


aE{L £ (T;x)} 


dt 


1 


+ “T P c,E{A,(T;*)} 

fTA(l/q) 


r^wq) f 

f Pi sup ^ £ ((1 - tg-sag) + ) > x, 

■/o Ug[o,i] 

6(1 - tq) V 6((! - tq- aq ) + ) < x j dt 


+ lim sup —— — —— 

40 aE{L £ (T;x)} 


r 

f p sup 6((1 - tq-saqy) > x, 

•lo Ug[o,i] 

6(1 - fg) V 6((1 -tq- aq) + ) < x j dt. 


Here lim sup of the left-hand side is finite, since the first term on the right- 
hand side is trivially finite, the second by (4.1) and the weak convergence 
used in (4.4), the third by (2.2) and (4.1) and the fourth and fifth, for a > 0 
small enough, by (2.1), (4.1) and (2.5). Hence the right-hand side limit in 
(4.6) is finite. Sending y { 0 and x | x in (4.7), and using (2.1), (4.1) and 
(2.5)-(2.6), we see that the following two limits exist and coincide: 

m > x} = o 

40 (q/q)E{L £ (T;x)} yio y 

In addition, it follows that k = sup TeI k(T) < 00 . Moreover, by inspection 
of (4.7), for Tel and x < x in J, 


(4.9) 


lim lim sup —Pi 
<40 £ | 0 q 


sup 4(1 ~tq)>x, 
ye[0,TA(l/g)] 

l((Tq)Al)/(aq)} ^ 

H {6(1 ~kaq) <x} l =0. 

k= 0 J 


By (2.7), (4.8) and (4.9), we have, for f G [0,1), x<x in J and Tel , 


lim inf lim inf — P 

40 40 h 


' [ii+ h )/q\ L(H'?)A1)/ (aq)\ 

U U {6 {n-kaq)>x} 

<n=\t/q\-\-l k=0 
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\h/n\ { L(( T 'i) A1 )/( a '?)J 1 

liminf liminf ---P< II {£ e (l — kaq ) > x} > 

MO MO h \ U 

Y l(.i+ h )/<l} L( i+h)/q][((Tq)/\l)/(aq)\ 

— lim sup lim sup — ^ ^ ^ 

£ J-° m= [£/ g j 4 -i n=m+l k,l =0 

P{C £ (m — kaq) > x, £ e (n — laq ) > x} 


(4.10) >k{T)F{x)+o( 1) 


— lim sup lim sup 

MO ejO 

gu ^ gP{g £ (l - ^g) > x, £, e (n - sq ) > x} 

r, sG [ 0 S TA(i/g)) P{£e(l - rq) > 0}P{£ e (l - sq) > 0} 

/L(( T 9)Al)/(a9)j P{ ^ e(1 _ kaq') > 0} \ 2 

*{ h «/9)P{&(l)>0 }) 

= k(T)F(x) + o(l) —> k(T)F(x) as a | 0, 


since (2.1), (2.2) and (2.5) show that the sum on the last row is bounded 
as e | 0. 

To finish the proof, notice that, given T E I and x £ J, (2.8) yields 


lim 

ej.0 


(4.11) 


r r LVfcj 

\[mh—h 2 )/q\ 

l((Tq)M)/{aq)\-l 

p n 

n 

n {&( n - ka q) ^ x } f 

f m=l n 

= L(m— l)h/q\+l 

k= 0 J 


-P 


■ l(h—h 2 )/q] L((Tg)Al)/(ag)J—1 

D fl {ie{n-kaq)<x} 

k =0 > 


[l/hh 


n =1 


= o. 


Using (4.11) together with (4.10), we readily obtain, for x E J, 


lim sup P < sup £ e (r) < x 


ej.0 


re[0,l/g] 


< inf lim sup 
~ re/ aj.0 


(4.12) 


lim sup ( 1 — h lim sup — P 
hlO y MO 


h 


■[(h-h 2 )/q\ L((Tq)Al)/(aq)J-l 

u u 

n=1 k =0 
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|i/M 


{P(n- kaq ) > x} 


< inf e - K(T)F(£) = e“ KF(i) -> e - Ki?(3:) as x j x. 

Tel 

On the other hand, (4.11) together with (4.8) and (4.9) similarly gives, for 
x G J, 


lirn inf P sup £ e (r)<x 
£ 1° Ue[o,i/ 9 ] 


> inf liminf 

TG/ cij,0 


(4.13) lim inf lim inf P < P) P) 

hlO elO 111 11 


\l/h\ [(mh-h 2 )/q\ [((T q) M) / (aq) ]-1 

n n n 

m=l n=[(m—l)h/q\+l fc=l 


d/ooj- 1 1 

Pi {PP - kaq) < x}> 

fc=0 J 


> inf liminf 

Tel ajO 


h. 


L((Tgf)Al)/(o9)J 


lim inf lim inf 1-P < 

/i|0 ej.0 y q 


(J {£,',s{n-kaq)>x}\\ 

k =0 J / 


Li/M 


> inf e - K(T)F(£) = e - Ki?(£) e ~ KF ^ as x T x. 

Tel 

Clearly, (4.12) and (4.13) add up to (4.1) or (4.2), depending on whether 
i = 1 or 2. □ 


5. Stable processes. We give two examples with strictly a-stable Levy 
motions, and prove (1.1) for totally skewed linear fractional a-stable motion 
(LFSM). 

Let {L(t)}t£«. be a separable a-stable Levy motion with L(t) spitivppp, 
where a G (0,2), /3 G [—1,1] and // G M, with g = 0 if a / 1 and (3 = 0 if 
a = 1. [See Samorodnitsky and Taqqu (1994) on basic facts and notation for 
a-stable processes.] 

It is tempting to try to extend Example 3 to a general SS Levy process 
L: 


Example 4. For (3^—1, Willekens (1987) showed that [see also Berman 
(1986)] 

(5.1) p( sup L{t) > nl ~ P{L(1) > -it} ~ bC a (l + (3)u~ a &su— >oo, 

Ue[o,i] J 
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where C a = (/ 0 °° x a sin(x) dx ) 1 . Hence, since L is SS with index l/a(l/a- 
SS), 


pj sup [L(t) - L q (t)} < e) 
Ue[o,i] ) 


= ( P<J sup L(t) < e 
^e[0,g] 


\ Li/gJ r l 

I P< sup L{t) < £ > 

) Ue[0,i-Li/?Jg] J 


~(l-P{L(l)>e/g 1 /“}) Ll/9j 

exp{ — ^C Q (1 + /3)e _ “} as g | 0 for e > 0! 

[It is an exercise in subexponentiality to use this in turn to recover (5.1).] 
Here L has jumps with positive probability. Hence the sampling problem 
does not really make sense: we are in fact computing the distribution of the 
largest positive jump. 

Further, again by (5.1), and in the case when a > 1 and /3 = — 1 also using 
that 

(5.2) p[ sup L(t) >u\ < for h, u > 0 

Ue[o,/i] J P{L(h)>0} 

[e.g., Albin (1993), equation 1.2], together with (5.3), we similarly get 

(l-P{|L(l)|>e/ (? 1 /a }) U/ 9 J 

>p( sup \L(t) — L q (t)\ < el 
Uelo.il J 


>(l-Pi sup L(t) >£/q 1/a \ -Pj inf L(t) < -£/q l/ ' 


te[o,i] 


1+1/9 


ie[o,i] 

(l-P{|L(l)|>e/g 1 /« }) U/ 9 J 

exp{— C a £~ a } as q | 0. 


For one case not covered by Example 4, that is, (1.1) for j3 = —l, Theorem 
1 applies: 


Example 5. Let a > 1 and /? = —!. Denote 


q(£) = 


where 


a/\ a 


2a In 


— (3 — 2a) In f 2 a In 


— 2 In I 


/ \/2ira \ 

\b* ,Xa y 


1 —a/(2\ a ) 


Aq; 

b 2 = _“_ 

“ 2(a — 1)| cos(7ra/2)| 2A “ 1 


and = 


a 


2(a — 1) 
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Recall that, as u —» oo [e.g., Samorodnitsky and Taqqu (1994), equation 

1 . 2 . 11 ], 


(5.3) 


P{S' a (fj,-1,0) >u} = P{5 q (1, —1,0) >u/a} 

P{N(0 ,b 2 a )>(u/a) x "} 


a 


Taking w = e/{2a\ a ln(l/e)) and q = q = aw/e (Q i = oo and Q 2 = 1), 
(5.3) gives (2.1) and (2.2) for L e = Li j£ , with F = G = e~'. J = M and I = 
( 0 ) 00). 

Let be a standard exponential random variable independent of L. Since 
{L e { 1 - qr)\L e {l - qr) > x) 4 £ + x, 
by (2.1), we now get (2.3) by observing that, as e j 0, 

{(Le ( 1 + q(t - r))\L e {l - qr) > x)} te(0 r] 
d f L(qqt) 




w 


+ (L e (l - qr)\L e (l - qr) > x) 


te(o,r] 


L 


cos 


/ na 


t) + x + C 


fe(o,r] 


Using (2.1) together with (5.2), we get (2.6) in the following way: 

P{sup r . e[TA( i /9 - )! (i /9 -) ) L e { 1 - rq) > x} 


lim sup ■ 

ejO 


P{L e (l)>0} 

P {Le{l ~ Tq) > x} 
eio' 1 ' P {L(q — Tqq) > 0}P{L £ (1) > 0} 

p -T-x 

0 as T —> 00 . 


< lim sup 


P{L(1)>0} 


Using 1/a-SS and (5.2), (5.3), (2.5) follows from the fact that, for y< 

y + fi, 


aP{L e 

< 


——-p( sup L e ((l - rq - sq) + ) > y + 6, L £ ((l - rq - aq) + ) < y\ 
(l)>bj Ue[o,al J 


—————T^P^ sup L e ((l-rq-sqy 
aP{L £ (l) > 0}^, Ue[o,a] 


- L e ((l - rq- aq ) + ) > 5 + k, 
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-M(l -rq~aq) + ) >y-k- l| 

< P{L(aqq ) > (<5 + fc)w}P{L g (l) >y-k- 1} 

“to «P{L(l)>0}P{L e (l)>0} 

< ^ P{L(|cos((7tq)/2)|) > {l/2){5 + k)/a l / a }F{L £ {l) > y - k - 1} 

~ to aP{L(l)>0}P{L e (l)>0} 

—» 0 uniformly in r 

as e [ 0 and a [ 0, since, by (5.3), 

P{L £ (1) > y - k - 1} 2k 

-P{L(1) > 0}- ~ K » e ■ 

Finally, (2.7) and (2.8) follow trivially from independence of increments. 

Now take a > 1 and (3 = —1, as in Example 5. Given an H £ (1/a, 1), 
consider the LFSM given by 

(5.4) {«()}.> 0 = {/ W + r)+) H - l l‘-(r+) H -V°\dUr)\ . 

UreR J t >0 

This process is H-SSSl and totally skewed to the left, £(t) ~ S a {a{f{l))t H , —1, 


Theorem 2. Let £ be given by (5.4). With 
s \ 1 / H 


1 (e) - (- 


1/Ao 


°m) 


H 


■ In 


HX n - 1 


H A 

— 2 In f 


In ( In 


V2 


ira 


\-I -l/(2ifA a ) 


and w(s) = He/(2\ a ln(l/e)) (c/. Example A), we have 

J s uPt e [o,i][£W -C(L*/9( £ )J9( £ ))] 


limP 

efO 


l w{e) 

= exp{— e~ x j for x € ’ 


< x 


Proof. With q{e) = q(s) = w{e) / {He), so that Q\ = oo and Q 2 = 1, 
(5.3) gives (2.1) and (2.2), for F = G = e - ', J = M and / = (0, 00 ), as in 
Example 5. 
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By Albin [(1998), Section 13], we have, with v(u) ~ cj.it 1 ^ a ^ andp(u) ~ 
C 2 v(u) 1 / h as u—* oo, for some constants ci,C 2 > 0, 


(5.5) 


P{£(1 ~p(u)t) > u+ (5 + x)u(u),£( 1) < u + xv(u)} 
< Ci<5~ d f e P{£(l) > u + xv(u)} 


for u large enough, t > 0 small enough and <5 > 0, for some constants C\, d > 0 
and e > 1. For any function p(u) = o(p(u )) as u —> oo, we thus have, by (2.1), 

P{£(1 -p(u)t) > u+ (5 + x)v(u),£(l) < u + xv(u)} 

hm - ——7 . ---= (J. 

u^oo P{£(1) >u} 

From this we get, by (2.1) and H- SS, with the change of variable u = e/q H , 
v(u) = w/q H ~ and p(u) = q/(l — q(r — t)) = o(u(u) 1// ^), for 

te(0,r), 

p {Ce(l ~ qr) > 5 + X, £ e (l - g(r - t)) < x} 

P{£e(l — qr)> x} 

c r+x 


P{C(1 )>e/q H } 


xpdi- 


qt 


> 


£ + w(S + x) 


(5.6) 


l — q(r — t)J q H [l — q(r — t)] H ’ 

£ + WX 


ea)< 


q H [ 1 — q(r — f)] 


H 


a r+x 


P{C(1 )>£/q H } 


xPdi- 


qt 


1 -q{r- 1 ) 


> 


£ + w(5 + x + r — t) 


7 if 


e(i)< 


£ + w(x + r — t) 


Q 


H 


r+a: P{£(l - pt) > u + (5 + x+ r - t)v,£(l) <u + (x + r — t)v} 


= e 


pm)> U } 


—> 0 as u-» oo. 

This in turn gives (2.3) with P r (t) = 1 since, by (2.1), 

liminf P{£ e (l — q(r — t)) > x|£ e (l — qr) > x} 
ej.0 


> lim inf lim inf 

<510 ej.0 


P{£e(l - q(r~t)) > x,£ e (l — qr) > 5 + x} 
P{Ce(l ~qr)>x} 
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> lim inf lim inf 

(510 ej.0 


P{g £ (l -qr)>5 + x} 

P{6(1 - &) > x} 


— lim sup lim sup 

(510 e|0 


P{6(i - q( r ~t)) < Xi£e(i~qr) > ^ + x} 
P{6(1 -qr) >^} 


= lim inf e 5 — 0 = 1 . 

510 

By (5.5) and H- SS, we get easily that [this is in essence what is going on 
in (5.6)] 

P{£(1 — pt)>u+(5 + x)v,£(l — ps ) 

< U + xv} < C2<5~ d (t — s) e P{£(l) > u + xv} 


for u large enough, t — s > 0 small enough, 5 > 0 and x and t in compacts, 
for some constants C 2 ,d > 0 and e > 1. Hence Albin [(1992), Proposition 2] 
(which does not really use stationarity) shows that, for x < x, uniformly for 
s > 0 in compacts, 


lim lim sup 

w—KX) 


PisuPtelo.a]^ 1 -p(s + t)) > u + xv,£(l - ps) < u + xv} 

«P{£(1)>«} 


= 0 . 


By the change of variable u = s/q H , and using H- SS, we now get (2.5), since 


lim lim sup 

aio £ ^o 


P{ su Pte[o,a]£e(i ~q{s + t)) >*,6(i - qs)<x} 
aP{6(l)>0} 


= 0 . 


Further, (2.6) follows easily from Albin [(1998), Theorem 12], together with 

( 2 . 1 ). 

By H- SS, with the new variable u [and since q(u) ~ C3/11], (2.8) follows if, 
given T > 0, 16 K and a,h> 0, for any si,..., Sk,t\,.. ., G [0,1/P{£(1) > 
it}] n U^Li(n — T/u,n] with Sj +1 — Sj > a/u, tj + 1 — tj >a/u and t\ — Sk > 

h/ PU(1)>«>, 


lim 

It—KX) 


{ k h' ' 

f]{£( s *) -£(L s d) <u + xv}. P|{^(ij) -C(L^'J) <u + xv} 

i= 1 3 =1 


(5.7) -p|f > |{C(si)-^(Ls i J)<u + xu}| 


x P { <u + xv} 


< 3 = 1 


= 0 . 


We show (5.7) by “truncating” the stochastic integral for £(f) — £(|_ij) from 
(5.4), 

{(i)-f(L‘J)= [ [((* + r) + ) f, - 1/ «_ ((LtJ+r) + ) «-i/« ]di(r) 

Jr£ R 








20 


J .M. P. ALBIN 


— + / 

J\r+t\<{l/2)h/V{Z(l)>u} J\r+t\>(l/2)h/P{£{l)>u} 

= h(t) + I 2 (t). 

Since k + k! < 1 + T/[aP{£(l) > u}], and by (2.1), it is enough to show 


that 

(5.8) 


limsupTAT^Uo 


for 5 > 0, 


u —>oc) aP{£(l) > u} 

because then we can subtract all I 2 (si) and I 2 {tj ) from the correspond¬ 
ing process values and, after having factorized the resulting independent 
probabilities, add the / 2 ’s again. The approximation error is asymptotically 
negligible, since by (2.1), 

lim lim sup (/c + A/)P{£(1) e [u + xv — 5v, u + xv + <5u]} = 0. 

^10 u —>00 

We get (5.8) from (5.3) using that (1 — x) H ~ l / a > 1 — x since, uniformly 
in t, 


cr{h{t)) a = j 


r+t|>(l/2)fc/P{£(l)>u} 


[(* + r) H - x ^ - ((LtJ + r)+)^- 1/a ] a dr 


/r+t>(l/2)/i/P{f(l)>u} 


(t + r) 


aH-l 


1 - 


1 - 


t - |tj 1 a 


t + r 


dr 


< / (t + ryn-^dr 

J r+t>{l/2)h/V{Z(l)>u} 

=o([p{((i)>«}rM). 

With the new variable u, using (2.1), (2.7) can be rewritten 
Lfc/P{£(i)>«}J 

lim lim sup sup Y F {t( 1 ~ r / u ) - £(i 1 ~ r / u \)> u i 

/ifO U^o o r,se[0,T) n= 2 


£(n — s/u) — £( [n — s/u \) > u] 


x [PU(l)>u}]- 1 = 0. 

By positivities present, this holds if 
LV p {£(i)>v}J 

limlimsup sup Y"' P{^(1)+^(n — s/u) — ^([n — s/uj) > 2u} 

hlO u >oo s g[o,T) n = 2 


x [P{e(l)>u}]" 1 =0. 
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When checking (5.9), we can disregard any finite number of terms, since by 
Minkowski’s inequality, and more or less immediate properties of LFSM, the 
scale 


<r(£(!)+ £("-«/“) -£(L n - s /“J)) 

< <r(£(l)) + a (£( n - s / u ) - £(l_ n - s/u\)) - 5 

= 2<r(£(l))-<5 

for r,sG [0, T] and n E {2,..., N}, for some 5 > 0. By (5.3), this ensures that 
\E p{{(l) +{(n - s/u) -{(|n-s/tij) > 2ii) 

h p {«i)>”} 

^. P{5„(2^(l))-^,-l,0)>2i.} 

- h p {fw > ”} 

—> 0 as u —> oo. 


Further, by routine estimates, the above scale is at most 

||((l + .)+)"- 1/ «_«.)+)"- 1/ « 

+ ((n-s/ii+-) + )" -1/ “ - (([n-s/“J +') + ) i,_ 1 / “ll„ 

< 2 1//a [<r(^(l)) + Kn H ~ l \ < 2 1 /“cr(£(l)) forse[0,T] and n > N, 


for u large enough and some constants K > 0 and Splitting the 

range of summation n £ {N, ..., |_h/P{£(l) > «}J} into two regions at n = 
[w/uj 1- ^, say, it follows readily, using (5.3), that also the remaining part 
of (5.9) to treat 


lim lim sup sup 

MO u->oo se[0,T) 


Lfe/P{e(i)>«}J 

E 


n=N 


P{£(1) +€(n - s/u ) - £(|_n - s/u\) > 2u} 


= 0 . 


□ 
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